TOWARDS THE KPP-PROBLEM AND logt-FRONT 
SHIFT FOR HIGHER-ORDER NONLINEAR PDES II. 
QUASILINEAR BI- AND TRI-HARMONIC EQUATIONS 



V.A. GALAKTIONOV 

Abstract. Extensions of ideas of Kolmogorov, Petrovskii, and Piskunov (1937) [35] on 
travelling wave propagation in the reaction-diffusion equation 

ut = Uxx + u{l — u) in M X IR+, uo{x) = H{—x) = {1 for a; < 0; for x > 0}, 

if (•) being the Heaviside function, are discussed. The present paper continues the study 
began in |23j for higher-order semilinear bi-harmonic and tri-harmonic equations 

Ut ^ -Uxxxx + u{l - u) and ut = u^xxxxx + u{l ~ u), etc. 

Here, some of the results are extended to the corresponding quasilinear degenerate 
parabolic models with nonlinearities of the porous medium type (n > 0), 

Ut^ -{\u\''u)xxxx+u{l-u), Ut^ {\u\'^u)xxxxxx+u{l-u), CtC. 

Two main questions to discuss are: 

(i) existence of travelling waves via any analytical/numerical methods, and 

(ii) their stability and derivation of the log i-shifting of moving fronts for some class 
of data uq (not for H{—x)). 



1. Introduction: the classic KPP-problem, its quasilinear extensions, 
and higher-order quasilinear parabolic pdes 

We begin with an introduction concerning classic results; see more details in [23l § 1]. 

1.1. The classic KPP— problem of 1937: convergence to TWs. In the KPP 
problem [38] (1937) 

(1.1) Ut = Uxx + u{l — u), X G M, t > 0; u{x,()) = Uq{x) in M, 
with the step (Heaviside) initial function 

(1.2) .„w = {J; III 
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the solution was proved to converge to the so-called mm«ma/ travelling wave (TW) solution 
corresponding to the minimally possible TW speed 

(1.3) the minimal speed of TW propagation : Aq = 2. 

Looking for a TW profile f{y) for arbitrary A > yields: 

u^:{x,t) = f{y), y = X — Xt, where / solves the ODE problem 

-A/' = r + /(I -f), ye R; /(-oo) = 1, /(+oo) = 0. 

This 2nd-order ODE, on the phase-plane {/, /'}, by setting /' = P{f), reduces its order: 



;i-4) 



(1.5) f = -A-/(^, 

and it was shown that there exists the minimal speed Aq = 2 and the unique (up to 
translation) minimal TW profile f{y)- Using the natural normalization 

(1-6) /(O) = i 

this minimal TW profile is defined uniquely. In addition, 

(1.7) /(?/)< in R. 

The characteristic equation for the linearized operator in (11. 4p has a multiple zero: 

(1.8) g" + 2g' + g = and g = e^'^ =^ (^ + 1)^ = =^ fii,2 = 
that yields the following asymptotic behaviour of f{y)'- 

(1.9) /(?/) = Coy e~^(l + o(l)) as y — )■ +oo, where Co > is a constant. 

Concerning the relation between the ODE TW problem (11. 5p for Ao = 2 and the PDE 
Cauchy one (11. ip . (II. 2p . the novel remarkable analysis in of convergence as t — )■ +oo of 
the solution of the Cauchy problem (II. ip . (II. 2p to the minimal TW (II. 4p was performed 
in the TW moving frame. This was very essential, and not in view of the obvious x- 
translational invariance of the equation (ll.ip : see below. Eventually, using PDE methods, 
the KPP-authors proved that the TW front moves like 

(1.10) Xfit) = 2t- g{t) ast^oo, with ^'(t) ^ 0, 
where the front location Xf{t) is uniquely determined from the equation 

(1.11) u{xf{t),t) = \ foralH>0. 
Then the convergence result of [3H] takes the form: 

(1.12) u{xf{t) + y,t) ^ f{y) as t ^ +oo uniformly in y G M. 
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1.2. Bramson's logt-front drift. The next important question is the actual behaviour 
of the TW shift g{t) in (11.101) for t ^ 1 (not addressed in [SB]). This is about an 
"centre subspace (manifold) drift" of general solutions of the KPP PDE (II. ip along a 
one-parameter family of exact ODE TWs {f{y + a), a € M}. 

This open problem was solved in 1983 by Bramson (see also [29]) by using pure 
probabilistic techniques (the Feynman-Kac integral formula together with sample path 
estimates for Brownian motion). It was proved that, within the PDE setting (II. ip . (II. 2p . 
there is an unbounded logt-shift of the moving TW front: 

(1.13) g{t) = k\ogt{l + o{l)) as t ^ +oo, with k = f , 

Thus, (I1.13P implies eventual, as t — )■ +oo, infinite retarding of the solution u{x,t) from 
the corresponding minimal TW (uniquely fixed by (II. 6p ). thought the convergence (I1.12p 
takes place in the TW frame. 

1.3. Known extensions: TW profiles for quasilinear second-order reaction- 
diffusion equations. The first attempts to extend the KPP results to quasilinear second- 
order react ion- diffusion equations 

(1.14) Ut = ikiu)u,,), + Q{u), Q(0) = Q(1) = 0, A;(n) > 0, 

includingthe porous medium diffusion operators with k{u) = m", n > 0, were known since 
the 1970l; see results and extra references inPEll21[I71[ISll3niEIlE21IMllMlEZlllIlll2], 
[H p. 34], [laiinillH], etc. 

However, it seems that any results concerning possible logt-front shift for quasilinear 
KPP problems including those with the degenerate operators of the PME type 

(1.15) Ut = iu^^%, + Q{u), Q(0) = g(l)=0, n>0 {u > 0) 

(with, say, positive non-Heaviside data; then + log t-shift may occur) were still unknown. 
Recall that, since (11.151) describes finite propagation, for the Heaviside initial data, the 
minimal travelling wave is indeed stable (just by a standard barrier-comparison approach, 
this was known since 1980s; see some references above). In this stability sense, the 
quasilinear KPP-2 problem (11.150 is simpler than the classic one for n = 0. However, 
we claim that logt-shift is possible for (11.150 for some classes of initial data. Indeed, 
by obvious reasons, in order to avoid a trivial comparison from below with a sufficiently 
shifted minimal TW f{x — Xot + a), with a 1, one needs that uo{x) must intersect all 
such TW profiles f{x + a) (evidently, the Heaviside one does not do that). Actually, this 
means that, for (11.151) . a "correct" KPP-type setting (a really difficult one) assumes that 
initial data do not have finite interface and are specially distributed, that do not allow to 
use straightforward upper/lower barriers for TW-like estimates above and below. 

""^It seems, the first ever sufficiently general existence TW result for (|1.14[) was obtained by A. Nepom- 
nyashchy, which was published in Russian in the technical Journal "Voprosy Atomnoi Nauki i Tehniki" 
( "Proceedings of Atomic Science and Technology" ) and was not translated into English. 
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If fact, it is easy to construct a number of explicit TW solutions for the PME with 
special source terms in (11.151) using the pressure variable: 



(1.16) = f =^ vt = {n + l)vvxx + {vxf + q{v), q{v) = nv"«^ Q{v«). 
The ODE for TW profiles f{y) then takes the form: 

(1.17) -Xf = (n + l)/r + ^ ifr + q{f), 

and, as usual, we are looking for a solution / having finite interface for y > and / — ?■ 1 
as y —7- — oo, with a linear spatial behaviour at the interfaces (this well corresponds to 
Darcy's law of their propagation for such weak solutions). 

Example 1. There exists the following explicit solution of (11.171) ((■)+ denotes the 
positive part): 

.^^g. m = ^, Ao = ^, where 

q{f ) = X{f - 1)2 - 2(n + 1)(/ - 1)=^ - (2(n + 1) + n±i)(/ - 1)1 

Example 2. There exists the following explicit solution of (I1.17P : 

/(?/) = ^^StI^, ^0 = ^, where g(/) = i (1 - /)[A(1 + /) 



'1.19) 

-{n + 1)/(1 - /) + (n + 1)/(1 + /)2 _ 1 H±i (1 _ /)(i + ff 



We claim that, in the corresponding PDE framework, such explicit TWs admit a ± logt- 
front shifting along an "affine centre subspace", for suitable classes of initial data (+logt 
for positive data, and — logt - for changing sign one). The derivations is similar (and 
simpler) to that in Section [9] for the quasilinear bi- and tri-harmonic flows. 

1.4. The present goal: extensions of the KPP—problem to higher-order quasi- 
linear parabolic PDEs. The main goal of the present research comprising [231 [21], and 

the present paper, is to show that the KPP-ideology can be extended to a variety of other 
more complicated higher-order semilinear and quasilinear PDEs with source-type terms. 

Thus, our goal here is to continue such a study, and, using various analytic, formal, and 
numerical methods), to show that such a general viewing of the KPP ideas make sense, 
and that many higher-order PDEs inherits some (but never all, of course) key features of 
this classic analysis. 

Let us note that, in [8[ [9] (see also a large list of references therein, as well as those in 
the present paper), TW profiles were studied to a class of quasilinear thm film equation^ 

(1-20) ut + {f{u))x = {h{u)ux)x - {c{u)uxxx)x- 

In particular, for more general TFE-type equations in X M , , stability of ID TWs was 
studied in [35] , though, as in most of the papers mentioned around, the authors dealt with 
nonnegative solutions of a free-boundary value problem, rather than the Cauchy one. 



^Indeed, this quasilinear model is most relevant to the present study of p.2ip and p.29p . though 
our solutions of the Cauchy problems are not nonnegative, are oscillatory, and of changing sign near the 
right-hand finite interface. 

4 



Thus, firstly, we will discuss some aspects of KPP-type problems for higher-order quasi- 
linear partial differential equations (PDEs), with the same Heaviside initial data. From 
applications, such KPP-type problems deal with nonlinear higher-order diffusion operators 
leading to well-known nowadays models; see references on various fourth and 2mth-order 
semilinear and quasilinear PDEs in[8l|9l[ni[Tl[T6l[l9l[20]. 

To this end, in the present paper, we extend some of the KPP-type results to the 
quasilinear bi-harmonic equation with the diffusion operators of the fourth-order porous 
medium type 

(1.21) Ut = —{\u\'^u)xxxx + u{l — u) in M x M_,_, where n>0. 

The corresponding TW with the speed of propagation A is then governed by the following 
fourth-order ODE: 

(1.22) u.{x,t) = f{y), y = x-Xt =^ _a/ = -(|/r/)/(4) + /(i _ /), 
with the singular boundary conditions at infinity: 

(1.23) f{y) =0 for ^ 1 and f{y) — )• 1 as y — )► — oo exponentially fast. 

The first condition in f ll.23p just takes into account that the degenerate equation fll.2ip 
describes processes with finite propagation of disturbances, a very well known fact for 
such PDEs; see key references and results in [22] . Moreover, solutions are oscillatory and 
changing sign close to finite interfaces. Therefore, the following example of a nonnegative 
TW solution cannot be generic: 

Example 3: nonnegative smooth TW for a bi-harmonic flow. We take n = 1 and 
consider the fourth-order parabolic equation with a TW / > 0: 

(1.24) ut = -{u').... + q{u) =^ -Xf' = -{fT^+q{f). 

For a special choice of the reaction q (see below), it admits the following explicit solution: 

(1-25) m = [^]\ 

where the non-oscillatory behaviour f{y) ~ (—y)^ close to the interface ?/ = 0~ is the 
actual decay for weak solutions of the Cauchy problem for n = 1; see Section HI It follows 
from fOSj) that 

(1-26) {-y) = j^-l, 

so that, substituting fll.25p into the ODE in fll.24p . performing all the differentiations and 
eventually expressing {—y) in terms of / via f ll.26p . one obtains the actual source term 
q{f), for which this is a solution. The required value of the speed A is then obtained from 
the asymptotic analysis near the interface, as y — )■ 0~, as follows: 

(1.27) f(y) = (-yf(l + Oiy)) =^ 3Xy' = -iy'p + qif ) + ..., so Ao = -120 < 0. 

Indeed, for other values of A, when the terms 0{y^) are not cancelled in the equation in 
fll.27p . one obtains a sufficiently "singular" reaction term q{f): 

(1.28) '?(/) ~ (-1/)' ~ as 2/^0", 
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i.e., q{f) is not Lipschitz continuous at / = 0, so that, in the PDE setting, there occurs 
a typical problem of non-uniqueness of solutions (if q{f) is a source term for / ^ 0). 

However, we will show in Section H] that the nonnegative TWs such as ( I1.25P are not 
generic in the sense that a. a. solutions are oscillatory and changing sign near finite 
interfaces. In other words, nonnegativity of solutions is not an invariant property of 
such quasilinear bi-harmonic flows in the Cauchy problem setting (but can be in a FBP 
one, what was first shown by Bernis-Friedman |6] for quasilinear thin film equations). 
Therefore, we must neglect any opportunity to deal with nonnegative solutions of such 
equations with arbitrary sufficiently smooth coefficients. 

For any A 7^ 0, the problem fll.22p . fll.23p is of the elliptic type, but it is not variational. 
Therefore, as in 121], one cannot use advanced methods for higher-order ODEs with 
potential operators associated with homotopy-hodograph and other approaches |13 | |36 | ITF] 
and/or Lusternik-Schnirel'man and fibering theory [261 127]. Therefore, the ODE f ll.22p . 
though looking rather simple, and, at least, simpler than most of related fourth-order 
ODEs already studied in detail, represents a serious challenge and cannot be tackled 
directly by known tools of modern nonlinear analysis and operator theory. 

Finally, we study TWs for the quasilinear tri-Harmonic equation 

(1.29) Ut = {\u\"'u)xxxxxx + u{l — u) in M x ]R+, where n > 0. 

1.5. Other related results of the present research. The previous paper [23] is de- 
voted to the semilinear parabolic equations for n = 0, i.e., 

(1.30) ut = -uxxxx + u{l-u) in MxM+ =^ _a/ = -/W + /(i _ /), 
and the boundary conditions now read 

(1.31) f{y) — > and f{y) — )■ 1 as ?/ — )■ —00 "maximally" exponentially fast. 

This "maximal" decay of f{y) at infinity somehow includes some kind of the remnants 
of a "minimality" of the possible TW profiles, though any direct specification of such a 
property is difficult to express rigorously. 

Thus, the present paper and [23] deal with higher-order parabolic equations of reaction- 
diffusion type, while, in the third part of this research [21], we study KPP-type problems 
for other types of PDEs including dispersion., hyperbolic, and other ones. Namely, in [21] , 
we will deal with higher-order hyperbolic and dispersion equations such as 

(1.32) utt = -Uxxxx + u{l - u) and Uut = D^x'^'^u + u{l - u) , etc. 

As for more "exotic" PDE models, as a formal but quite illustrative example, we consider, 
in [21], higher-order dispersion equations and end up with the following one: 

(1.33) Dlu = D]^u + u{l-u), 
with eleventh-order ODE for the TW profiles 

(1.34) - aV(9) = /(II) + /(I - /) in M (plus ([I23D). 

We also consider in [21] an example of a quasilinear dispersion equation with a similar 
nonlinearity. 
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1.6. Main questions to study. Thus, for quasilinear KPP-problems fll.2ip and fll.29p . 
the main questions to study, here and in [23], [21] (though in the quasilinear cases the 
results achieved are weaker), are: 

(I) The problem of TW existence: existence of travelling waves via analyti- 
cal/numerical methods. Firstly, we easily show the positivity: A > always. 

(II) The logt-SHlFT problem: given some speed Aq > 0, derivation of a possible 
log t-shifting of the moving front in the problem (11. ip . for classes of initial data uo{x) (for 
Heaviside data (11 ■2p . we expect no shift and convergence to a TW; cf. some convincing 
stability results in [H [9l [35] for thin film equations). As in the semilinear cases, this 
phenomenon is also connected with a kind of an "(affine) centre subspace behaviour" for 
the rescaled equation. 

We thus show that log t-shifting phenomena, for some classes of data) are also available 
for quasilinear degenerate KPP-4 and other problems, though is less generic than for 
semilinear higher-order equations, since require, in general, still unknown hypotheses on 
Uo{x). 

Unlike the semilinear analytic PDEs for tt, = in [23] , the rescaled quasilinear equations 
for ?2 > are not analytic, so we cannot proceed with a deeper study of the omega limits 
for such equations. 

2. The basic higher-order KPP-4n problem 

Consider the basic KPP-(4,l)n (or simply KPP-4n, that cannot confuse in the para- 
bolic case) problem (ll.30p . and let us begin with its ODE counterpart (ll.22p . (I1.23p . 

2.1. A is always positive. As in [23| § 2], we first prove a first simple result on the 
positivity of admissible speeds A: 

Proposition 2.1. 

(2.1) // there exists a solution f{y) ^ of ( ^22^ . ( K^^ . then A > 0. 

Proof. Multiplying the ODE fll.22p by d/l"/)' and integrating by parts over M and using 
the boundary conditions (I1.23P yield 

, , -Kn + l)I\fnfr=[F{f{y))]t^ = -F{l)<0, 
2.2 , 

where F{f) = (n + 1) |z|"z(l - z) dz. □ 

2.2. Numerical construction of TW profiles. Again, as in [23], we begin with pre- 
senting first numerical results, which directly show the global structure of such TW profiles 
to be, at least partially, justified analytically. We use the bvp4c solver of the MatLab with 
a sufficient accuracy; see more details in [23], § 2]. Note that, as the initial data for further 
iterations, we always took the Heaviside function 

(2.3) initial data for numerical iterations are often H{—y), 
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Figure 1. The TW profiles F{y) of ([MD, (HTZ^ for A = 1, n = 0, 0.1, 0.2. 

i.e., as in (11. 2p . This once more had to help us to converge to a proper "minimal" profile 
(indeed, there are many other TW profiles), though, of course ,this was not guaranteed a 
priori. We kept this rule for all other KPP-(A;, /) problems of interest in [231 121] ■ 
For convenience, we perform all the calculations for the function 

(2.4) F{y) = \f{yTf{y) =^ F(^^ = ^ iFp^F' + |Fr^F(l - iFp^F), 

so that the ODE in (ll.22p becomes semilinear. 

Figure [T] shows TW profiles F for A = 1 for sufficiently small n = 0,0.1,0.2. The 
semilinear case t?, = was used to get a good comparison with the results in |23l § 2] . For 
larger n = 0.8, with the same A = 1, the numerical results are shown in Figure [2j 

We next fix n = 1 and A = 0.5 as Figure [3] shows. A more detailed structure of 
oscillations of solutions about equilibria F = (a) and F = 1 (b) is presented in Figure 
m Separately, Figure |5] shows TW profiles F{y) as in (12. 4p for n = 1 and A = |. 

Finally, we must admit that, for n > 0, and, especially, for larger values n ~ 1, it is much 
more difficult to get reliable well- converging numerical results unlike the more straight- 
forward semilinear case n = 0, [23]. Many our attempts, even rather time-consuming 
(lasting for several hours of the MatLab) led to oscillatory non-converging profiles and/or 
to singular Jacobians (a kind of "blow-up" of computations). 

However, our numerical experiments produced a rather convincing evidence of existence 
of TW profiles for a variety of > and velocities A > 0, which we will need to develop 
further. 
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Figure 2. The TW profiles F{y) of ([231), for A = 1, n = 0.8. 
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Figure 3. The TW profiles F{y) of dH, fOSjl for A = 0.5, n = 1. 
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KPP-4n: -X f= - (tfri)'*' + f(1-f), TWs for 1=0.5 and n=1 




(a) oscillations about F = 

F(y) KPP-4n: -X f'= - (Ifrt)'"' + f(1 -f), TW profiles for 1=0.5 and n=1 
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(b) oscillations about F = 1 

Figure 4. An oscillatory convergence of the TW profile F{y) to 0, y ^ 1 (a) 
and to 1, y <C —1; A = 0.5 and n = 1. 

Thus, we begin with necessary local analysis of the behaviour of TW profiles f{y) close 
to equilibria / = 1 and / = (see the next section), in order to fix dimensions of their 
stable and unstable manifolds, for further matching. 

3. The 2D stable bundle as y -oo and instabilities 
This analysis is not much different from that for n = in [221 § 2.3]. 



3.1. Linearization about / = 1. Thus, setting f = 1 + g in fll.22p and linearizing yield 
the following characteristic equation: 

(3.1) -Xg' = {n+l)g^^^+g and g{y) = e^^ =^ H_{fi, A) = (n+l)//^-A//+l = 0. 
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Figure 5. The TW profiles F{y) of ([23D, for A = 0.25, n = 1. 

Therefore, for A = 0, we have 2D stable (Re(-) > 0) and unstable (Re(-) < 0) bundles 
with the roots 

(3.2) ^^(0) = ±^(r2 + l)-i and /i±(0) = ±i^(ri + 1)^3. 

By continuity, for all small A > 0, there exists 2D stable manifold of the equilibrium 1 
with the roots 

(3.3) stable bundle : /^+(A) and /i+(A). 

Proposition 3.1. For any n > 0, at least for all small |A| > 0, the linearized equation 
(13. ip {and hence the KPP-An one in (ll.22p for / ~ 1 for y ^ — 1) admits a 2D sta- 
ble family of oscillatory solutions as y —oo, and a 2D unstable one of exponentially 
divergent orbits. 

3.2. Local blow-up and other instabilities to — oo. In order to verify the global 
continuation properties of stable bundles, one need to check whether the nonlinear ODE 
( I1.22P admits blow-up and the dimension of such an unstable manifold. To this end, we 
re-write it down and, as usual, neglect the linear lower-order terms, which by standard 
local interior regularity are negligible for |/| 3> 1: 

(3.4) (|/|"/)(4)=_/2^^^^_^.^__^2(i^^(^)) f^^_ 
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For any n G [0,1), the unperturbed equation has the following exact blow-up solution: 
for the function F in fl2.4p . 

(3.5) fW = -|F|^ =^ Fo{y) = ~Coin){Yo-yr^ ^-oo as y ^ Y,' , 

where Iq ^ 1^ is a fixed arbitrary blow-up point and Co{n) > is an easily computed 
constant. Studying the dimension of its stable manifold, as in [231 § 2.4], yields 

Proposition 3.2. For the ODE (11.221) . with any n G [0, 1), the stable manifold of blow-up 
solutions is ID depending on a single parameter being their blow-up point Yq G M. 

For n = 1, the asymptotic ODE (13. 5 p has the form 

(3.6) F(^) = -|F| =^ Fo{y) = -e'^ ^ -oo as y^ -oo, 

so, instead of having finite y-blow-up, one gets this negative exponential growth for y <C 
— 1, with a similar one dimension. 

Finally, for n > 1, we have algebraically growing solutions: 

(3.7) Fo{y) = -C,{-y)^ ^ -oo as y ^ -oo (Q > 0). 

4. Oscillatory solutions near finite interfaces: 3D asymptotic bundle 

This part of the asymptotic analysis is essentially different from that for the semilinear 
case n = in [231 § 2]; however, such results are already known [22j, so we can omit some 
details. 

Thus, for n > 0, we describe the generic oscillatory behaviour of solutions of (12.41) close 
to finite interfaces. According to pioneering results in [5l [7], since 1988, it was known 
that ODEs like (12. 4p admit compactly supported solutions of changing sign near finite 
interfaces. Let f{y) vanish at the interface y ^ yo > 0, so that f{y) = for y > yQ. 
Then, making for convenience the reflection y ^-^ yQ — y, with y > small enough, and 
keeping the leading first two terms in (12.40 for y ^ 0^, after integration once we obtain 
an exponentially small perturbation of the following third-order equation (the sign "— " 
on the right-hand side appears because of the reflection in y): 

(4.1) F'" = -AlFp^F, y > 0, F(0) = 0. 
We next scale out the positive constant A to get the ODE 

(4.2) F'" = -\f\ y>0. 

We next describe oscillatory solution of changing sign of the ODE (14. 2p . with zeros 
concentrating at the given interface point y = 0~^. Let us mention again that oscillatory 
properties of solutions are a common feature of higher-order degenerate ODEs. We refer 
to first results in [5l [TJ |22], to [TU [15] (thin film equations), and to [2H1 Ch. 3-5], where 
further examples can be found. 

To this end, by the scaling invariance of (14. 2p . we look for its solutions of the form 

(4.3) F{y) = y^v{s)-, s = \ny, where = -^^^^^-^ > 3 for n > 0, 
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where (p{s) is the so-called oscillatory component. Substituting ( I4.3p into (14. 2p yields the 
following third-order equation for (^{s): 

(4.4) P-siif) = -\ip\-^^ip, 

where Pk denote linear differential polynomials obtained by a simple recursion procedure 
(see [281 P- 140]), so that 

P,{ip) = ^' + fi^, P^iip) = ^" + (2/i - l)ip' + /i(^ - 1)^, 

P^iip) = if'" + 3(^ - l)ip" + (3^2 _ + 2)ip' + 12(^1 - l)(/i - 2)ip. 

According to (14. 3p . we are interested in uniformly bounded global solutions (p{s) that 
are well defined as s = Iny — )■ — oo, i.e., as y — > 0+. The best candidates for such global 
orbits of (14. 4p are periodic solutions (p*{s) that are defined for all s G M. Indeed, they can 
describe suitable (and, possibly, generic) connections with the interface at s = — oo. See 
[22] for the following and other related results. 

Proposition 4.1. For any n > 0, (14.41) has a periodic solution <y9*(s) of changing sign. 

Two problems remain open: 

(i) uniqueness of the periodic solution (p*{s), and 

(ii) stability (hyperbolicity) of V9*(s) as s — )■ +oo. 

Numerically, we have obtained positive answers to both questions. In particular, (i) and 
(ii) imply that there exists a unique (up to translation) periodic bounded connection with 
s = — oo, where the interface is situated. 

The convergence to the unique stable periodic solution of (14. 4p is shown in Figure [6] 
for various n > 0. Different curves therein correspond to different Cauchy data if{0), 
if'{0), f"{0) prescribed at s = 0. For n smaller than |, the oscillatory component gets 
extremely small, so an extra scaling is necessary, which is explained in [14, § 7.3]. A more 
accurate passage to the limit n — )■ in (14. 4p is done there in Section 7.6 and in Appendix 
B. This explains the continuous deformation as n — )■ of oscillatory structures in (14. 3 p to 
linear ones in the exponential tail of the kernel F{y) of the fundamental solution of the 
bi-harmonic operator Dt + -D^. 

In (d), we also present the periodic solution for n = +oo, where (14. 4p takes a simpler 
form (see an algebraic construction of the unique periodic solution in ^4, § 7.4]) 

P-iiv) = -signv?. 

Finally, given the periodic f^:{s) of (14. 4p . as a natural way to approach the interface 
point yo = according to (14.30 . we have that the ODE (14. 2 p generates at the singularity 
set {/ = 0} 

(4.5) a 3D local asymptotic bundle with parameters yo, a phase shift in s i— )■ s + Sq, 
and the parameter £ > of the scaling group for the ODE (14. 2p . 

3(n + l) 

(4.6) F^e^^F, y^ey (5 > 0). 

Notice that this scaling invariance has been lost in the approximate ODE (14. 2p . 
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(a) n = 0.75 



(b) n = 1 




Remark: nonnegative solutions for A < (Example 3). One can see that f l4.4p 
does not admit nonnegative solutions and equilibria. However, if A < in (14. ip . then 
scaling out such A leads to ODEs fl4.2p and f l4.4p with the opposite sign on the right-hand 
side, so that there exists a positive equilibrium 

(4.7) PM = \^\-^^cp =^ ip,= [^(^ - 1)(^ - 2)]-^ 

(the rest of solutions are still oscillatory). Existence of such a particular positive solution 
in (14. 7p is the "local" origin of the global nonnegative TW in Example 3 in Section [H 
though such a behaviour is not structurally stable for both ODEs involved. Note also 
that the calculation (12. 2 p does not apply to (ll.24p with a more complicated reaction g(/). 



4.1. Well-defined matching of stable manifolds of equilibria and 1. Similar to 
[23] for n = (when the conclusion is indeed more straightforward), we thus observe that, 
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for n > and sufficiently small A > 0, there is a well defined matching procedure of stable 
manifolds of equilibria / = (see (14.51) ) and / = 1 (Proposition 13. 11) . The overall relation 
of the dimensions is 

(4.8) 3(j^»i,/^o) - 1 = 2(j;«_i,/^i), 

where "—1" stands for the dimension of the unstable (blow-up for n < 1) manifold in 
Proposition 13.21 In other words, (14. 8 p implies a well posed algebraic system of two equa- 
tions with two unknowns. In the case of analytic manifolds as in [23], this guarantees at 
most a countable number of solutions (or a finite one of uniformly bounded TW profiles). 

In the present degenerate case n > 0, any analytic dependence on parameters is also 
plausible but difficult to prove. However, we think that (14. 8 p conffims that, for small 
A > 0, the family of TW proffies is always discrete (what have been seen in numerical 
experiments), and then such (most probably, finite) number of A-branches can be extended 
for larger A > by classic nonlinear integral operator theory [211 [12] • Note that a fully 
global extension in A is hardly possible: as was shown in [23] for n = 0, there is a clear 
phenomenon of existence of a maximal speed Amax such that, for A > Amax, TW proffies 
f{y) are nonexistent. By a continuity argument, we expect that a similar phenomenon 
takes place for n > 0, but this is difficult to check even numerically. 

5. A FEW WORDS ON "MORE QUASILINEAR" KPP-4?7, PROBLEM 

This is about the KPP-setting for the following quasilinear equation: 

(5.1) = -{\uruu,,+\uru{i-\uru) =^ -a/ = -(|/|"/)W+|/r/(i-i/r/), 

where the source contains the same nonlinear term \u\"'u as the porous medium diffusion 
one. A one advantage of this model is that the corresponding ODE for F = I/I"/ is less 
complicated than that in (12. 4p : 

(5.2) FW = _A_|F|-;rrTF' + F(l-F), 

which, as might be expected, could improve convergence of numerical schemes. However, 
this did not happen, and convergence, though being slightly better, was not improved 
somehow essentially. 

In Figures [7] and IH], we show TW proffies F{y) for n = 0, 0.5, and 1 for the two cases 
A = 1 and A = 0.5. 

An unusual application of such a quasilinear bi-harmonic equation (15. ip is that it allows 
to treat the "fast diffusion" case n < 0; more precisely, n G (—1, 0). In Figure we show 
such oscillatory profiles for n = 0, —0.25, and —0.5 for A = 0.5. The convergence for n < 
is much better than in the standard porous medium case n > 0. It is worth mentioning 
that, since for u ^ the source term becomes a non-Lipschitz function, Q{u) ~ +|m|"'m, 
n < 0, there appears a standard non- uniqueness of solutions even in the simple ODE: 

(5.3) ut = \u\''u, t>0, u{0) = =^ 3 u{t) = {\n\t)\^ > 0, 

together with the trivial solution u{t) = 0. In such cases, in PME and fast diffusion theory 
with blow-up, extinction, quenching, and other singularities, an extended semigroup of 
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F(y) KPP-4n: -Xf'= -(|fr f)<'*'+|f|''f(1-|f|"f), TW profiles for X=1 , n=0, 0.5, 1 
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Figure 7. The TW profiles F{y) of (Ol for A = 1, n = 0, 0.5 and 1. 

'"(y) O'^^'+lfl^fll-lfrO. TW profiles for;L=0.5, n=0, 0.5, 1 




Figure 8. The TW profiles F{y) of (jEZD, (Ol for A = 0.5, n = 0, 0.5, and 1. 
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F(y)KPP-4n: -X f'= -(Ifl" f)<'*'+|frf(1-|frf), TW profiles for;i=0.5, n=0, -0.25, -0.5 
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Figure 9. The TW profiles F{y) of (IE2]), fOHl) in the "fast diffusion" 
case for A = 0.5, n = 0, —0.25 and —0.5. 

proper extremal (maximal or minimal) solutions is constructed; see [2T[ Ch. 6,7] and 
references therein. However, such a theory is fully nonexistent for higher-order parabolic 
flows with such singularities. Creating such a theory for equations with no Maximum 
Principle is a hard open problem. 

6. A FEW WORDS ON NON-DIVERGENT KPP THIN FILM PROBLEM 

We next briefly review some results for the KPP problem for the parabolic thin film 
equation (TFE-4) (n > 0) 

(6.1) ut = -(kr«...). + nil -u) =^ -xf = -d/rn' + /(i - /)• 

In Figure [TOl we show TW proflles for A = 1 and sufficiently small n = 0.3, 0.6, and 0.9, 
and, for comparison, we present the "linear" one for n = (i.e., for the KPP-4 problem 
in [23] )• Next Figure [TT] shows profiles for larger n = 1.5, 2, 3, and 4. 

In Figure [121 we show an enlarged image of the behaviour of f{y) from Figure [TT] for 
large y G [10,17]. It is seen that f{y) becomes less oscillatory close to interfaces as n 
reaches about 2. We believe that a changing sign feature for n > 2 in this figure can be 
related to the necessary regularization of the degenerate thin film operator in (16.11) . where 
we replace 

(6.2) |/|"^(e2 + /2)t, with e = 10-^ 
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Figure 10. The TW profiles f{y) of fOSjl for A = 1, n = 0, 0.3, 0.6, 
and 0.9. 

since smaller e = 10"'^, etc., always led to divergence of the numerical scheme; see below. 

Note that the convergence for larger n got very slow; for n = 4.3, we still got some 
profile, which is close to that for n = 4 in Figure [TTl while, for n > 4.4, "a singular 
Jacobian" (a full non-convergence) always appeared. 

The oscillatory behaviour of f{y) near the interface was studied in [TU [25]. It was 
shown that a periodic oscillatory component ip{s) (a full analogy of that in (14. 3p . ( 14. 4p ) 
exists up to a critical homoclinic bifurcation exponent nh, 

(6.3) < n < nh e (|, n+), where n+ = = 1.9019238... . 
Numerically, rih is given by 

(6.4) rih = 1.75987... . 

Therefore, it is expected that, for n > nh, f{y) exhibits either a finite number of zeros 
close to the interface, or even becomes nonnegative (at least, for n >2). 

Finally, existence of a logt-shift of moving front for the PDF KPP-TFF-4 problem 
(16. ip is shown precisely in the same way as in Section |9l 
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Figure 11. The TW profiles f{y) of fOSjl for A = 1, n = 1.5, 2, 3, 
and 4. 

7. The quasilinear KPP-(6,l)n problem 

We now, more briefly than above, consider the KPP-(6,l)n problem (I1.29P and its ODE 
counterpart 

(7.1) -A/ = (i/r/)(^) + /(i-/), F = i/r/, 

with singular boundary conditions ( ll.23p . 

Figure [13] shows TW profiles for A = 1 and n = 0.25. For comparison, we also put the 
"semilinear" profile for n = 0. This confirms a rather non- surprising fact that there exists 
a continuous dependence of f{y) on n — )• 0+ (this even can be proved for such ODEs). 

In Figure [TU we present a similar comparison for the more oscillatory case A = 0.2, 
where we again observe existence of a clear "homotopy" deformation as n — > 0. 

A proper dimensional analysis of the linearized bundle as y — )■ — oo (i.e., as / — t- 1) is 
performed similarly, as in |23] . 

The oscillatory periodic-like behaviour at the finite interface as y (i.e., as / — )■ 0) 
is performed as in [15]; see also [281 P- 142]. 
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f (y) X 1 0"'* KPP-TFE-4: - X f=-{|fr f"')'+f(1 -f). . n=0, 1 .5, 2, 3, 4 
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Figure 12. Enlarged behaviour of TW profiles f{y) close to zero of (16. ip . 
^L23\f from Figure [IH A = 1, n = 1.5, 2, 3, and 4. 

Note that the corresponding blow-up problem for n G (0, 1) (cf. the ODE (13. 5p ) 

(7.2) (i/r/)(^) = f 

is now more difficult since (17. 2p admits oscillatory solutions, which can be studied as in 
[m § 7] and in [15] by introducing an oscillatory component represented by periodic or 
other functions. 

8. Very briefly on KPP-8n 
In the KPP-St?,, we deal with the following ODE problem: 

(8.1) Ut = -Dl{\uru) + u{l-u) =^ (ODE) - A/' = -(|/r/)^'^ + /(I - /)• 

Figure fT5l shows TW profiles F for A = 0.5 and n = 0.5, 1, 2, plus, for comparison, for 
n = from [23, § 4]. This shows a principal positive answer on the TW existence question. 
It is seen that oscillations about / = decrease as n increases, though finite interfaces 
for n > (unlike n = with an exponential decay as y ^ +oo) are obviously invisible. 

Stable bundles as / — )■ 1 (y — )■ — oo) are studied as in [23], while, close to finite interfaces, 
as / — )• 0, periodic oscillatory components for seventh-order ODEs are constructed as in 
[281 P- 142]. 
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Figure 13. The TW profiles F{y) satisfying (El), f Oal) for A 
n = 0.25, n = 0. 



1, and 



9. The origin of logt-SHiPT: centre subspace balancing 

We now explain the origin of possible log t-shiftings (say, retarding) from the TW for 
the PDE higher-order quasilinear parabolic KPP-problems (11.211) or (ll.29p . Note that 
the actual proof of such log t-drift assumes a delicate matching of the solution behaviour 
on compact subsets in the TW variable, i.e., in the Inner Region, with a remote Outer 
Region for y ^ 1, where the infiuence of the nonlinear term — is negligible, and the 
actual behaviour is governed by the quasilinear bi- (or tri-) harmonic operator. Such 
a matching, eventually, is supposed to describe classes of initial data, for which such 
(±) log t- shifting (or no shifting at all, say, a pure convergence to a TW in its moving 
frame). We do not present here this kind of a procedure of matching of asymptotic 
expansions of those two regions, which is extremely difficult in the quasilinear /degenerate 
case. Therefore, we restrict to an Inner expansion analysis. 

9.1. Linearization and rescaled equation. Thus, we consider a KPP-type problem 
for a quasilinear PDE 

(9.1) ut = A{u) + u{l - u) 



m 



X 



^+5 



with some proper step-like data uo{x), with, however, some positive, negative, or even 
oscillatory "tails" as y — )■ -|-oo, in order to create a necessary ±logt-shift. Here, in 
(19. ip . A(w) is a proper homogeneous isotropic quasilinear differential operator satisfying 
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KPP-6n: -X f= (|f|" f)<^'+f{1-f), TW profiles for>.=0.2, n=0, 0.25 




Figure 14. The TW profiles F{y) satisfying (El]), fOal) for A = 0.2, and 
n = 0.25, n = 0. 



some extra conditions specified below. We fix, as key examples, the quasilinear bi- or 
tri-Laplacian operators of the porous medium type 

(9.2) A(m) = or A{u) = Dl{\u\''u), n > 0. 
We assume that, for some fixed Aq > 0, the corresponding ODE problem 

(9.3) -Ao/ = A(/) + /(!-/), 

with the conditions f ll.23p admits a unique solution /. 

Attaching the solution u{x, t) to the front moving and setting, for convenience, Xf{t) = 
XqI — g{t), the PDE reads 

{9A) u{x,t) = v{y,t), y = x-Xot + g{t) =^ vt = A{v) +v{l -v) + XoVy - g'{t)vy. 
We next linearize (19. 4p by setting 

(9.5) v{y,t) = fiy) + w{y,t), 

that yields the following perturbed equation: 

Wt =Bw — g' (t) f — g' {t)wy — w'^ , where 

Bw = A'{f)w + {l-2f)w + XoWy and A'(/)w; = (n + l)A(|/|"w;). 

Assuming that, in this g{t)-movmg frame, there exists the convergence as in (I1.12p . so 
that w{t) — )■ as t — )■ +00, one can see that the leading non-autonomous perturbation in 
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KPP-8n: -Xr= f)<^'+f(1-f), TW profiles for;i=0.5, n=0, 0.5, 1 , 2 
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Figure 15. The TW profiles f{y) satisfying (EI]), f Oal) for A = 0.5 and 
n = 0, 0.5,1, 2. 

(19. 6p is the second term on the right-hand side. However, as we show, the last two terms, 
though negligible, will define a proper logt-shift of the front. 

9.2. Where logt-shift comes from. Again, as in the semilinear case n = 0, we note 
that the rescaled equation (19.61) is essentially non- autonomous in time, so we cannot use 
powerful tools of nonlinear semigroup theory; see |l0] . However, using a formal asymptotic 
approach, we will trace out some definite centre subspace behaviour after an extra rescaling 
and balancing of non-autonomous perturbations. 

Thus, as usual (see Introduction), we assume that g'(t) — )■ as t — )■ -|-oo sufficiently 
fast, i.e., at least algebraically, so that 



Under the hypothesis (19.71) . the only possible way to balance all the terms therein (in- 
cluding the quadratic one —w"^) for t 1 is to assume the asymptotic separation of 
variables: 

(9.8) w{y,t) = g'{t)ij{y) + e{t)(p{y) + where \e{t)\ <t: \g'{t)\ as t ^ +oo. 
Here, we omit higher-order perturbations. Substituting (19. 8p into (19.81) yields 



(9.7) 



/(t)|«|^'(t)| for t»l. 



(9.9) 



g"{t)ij + e'{t)v + ...= g'{t){B^-f) 

+ - {9\t))\i^' + - 9\t)e{t){v' + 2^^) - e\t)v' + ... . 
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Using (19. 7p and (19 .Sp in balancing first tlie leading terms of the order 0{g'{t)) yields 
the elliptic equation for ip: 

(9.10) 0{g'{t)) { = 0i])): Bij-f' = 0. 

Then balancing the rest of the terms in (19.90 requires their asymptotic equivalence, 

(9.11) g"it) ~ -ig'it)f ~ sit), i.e., git) = klogt, g'{t) = f , g"{t) = -A, e{t) = ^. 
Then, we obtain the second inhomogeneous singular Sturm-Liouville problem for ip: 

(9.12) 0{-^): Bip = k^lj + k'^{%l)' + '4)'^). 

Thus, the first simple asymptotic ODE in (19. lip gives the log t-dependence as in (ll.lSp . 
Finally, we arrive at the following system for {-0, y)}: 

= /', 



(9.13) 



Solving this system, with typical boundary conditions as in (ll.23p . allows then continue 
the expansion of the solutions of (19.60 close to an "affine (i.e., shifted via /' on the RHS) 
centre subspace" of B governed by the spectral pair obtained by translation in (19.30 : 

(9.14) Ao = and i;,{y) = f\y). 
The asymptotic expansion for t ^ 1 then takes the form 

(9.15) w{y,t) = \i,{y) + j,^{y) + ..., 

which can be easily extended by introducing further terms, with similar inhomogeneous 
Sturm-Liouville problems for the expansion coefficients. 

As for 72, = 0, B does not have a discrete spectrum, so we cannot get a simple algebraic 
equation for k by demanding the standard orthogonality of the right-hand side in the 
second equation in (I9.13P to the adjoint eigenvector ip^ of B* in some "weighted" metric 
(■, ■)* (in which the adjoint operator B* is obtained, if any), like 

(9.16) k: {ki{j + e{ij' + ij^),tl;*), = 0. 

Therefore, it seems, the system (I9.13P cannot itself determine the actual value of k therein. 
As we have mentioned, the latter requires a difficult matching analysis in Inner and Outer 
Regions, which, for the KPP-4n (and all other problems) remains an open problem. 
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